Instability of the rhodium magnetic moment as origin of the metamagnetic phase 

transition in a-FeRh 
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Based on ah initio total energy calculations we show that two magnetic states of rhodium atoms 
together with competing ferromagnetic and antiferromagnetic exchange interactions are responsible 
for a temperature induced metamagnetic phase transition, which experimentally is observed for 
stoichiometric Q-FeRh. Taking into account the results of previous and newly performed first- 
principles calculations we present a spin-based model, which allows to reproduce this first-order 
metamagnetic transition by means of Monte Carlo simulations. Further inclusion of spacial variation 
of exchange parameters leads to a realistic description of the experimental magneto-volume efi'ects 
in Q-FeRh. 
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I. INTRODUCTION 



In 1938 Falloti^ discovered that ordered bee FesoRhso 
undergoes a first order metamagnetic transition from an 
antiferromagnetic (AF) ground state to a ferromagnetic 
(FM) phase with increasing temperature. This transition 
occurs at w 320 K and is accompanied by a volume 
increase of about 1 The Curie temperature Tq of 

the FM phase is of the order of 670 K."'"'^ In contradic- 
tion to the first hypothesis of Fallot and Horcartf^ X-ray 
diffraction measurements showed that the transition is 
isostructuraliiSi From Mossbauer and neutron diffraction 
measurements one knows that the FM phase has coUinear 
magnetic moments of 3.2 hb per Fe atom and 0.9 fiB per 
Rh atomi At low temperatures an AF-II spin structure 
is found with Fe moments of 3.3 /iB and with vanishing 
Rh moments (see Figure^. Application of hydrostatic 
pressure suppresses the FM phase, i. e. for a critical pres- 
sure of 60 kbar the system immediately transforms from 
the AF to the paramagnetic (PM) phase. 

An early explanation for this behavior was based on the 
phenomenological exchange inversion modeM by Kit- 
tel (which originally was designed to explain metam- 
agnetic transitions in other materials like Cr-modificd 
Mn2Sb with a layered magnetic structure). In this 
model the exchange parameter varies linearly with the 
lattice constant and changes sign for a critical value a^. 
However, experimental findings like the rather large en- 
tropy chang8^i^*ii*i2ii^ at Tm, which is of the order of 
AS'(rm) w 12.5- 19.7 Jkg-^K-i, as well as elastic prop- 
erties could not correctly be describedii^ Tu et alm^ 
used a different approach by considering the large dif- 
ference of the low-temperature specific heat constants, 
7, of the AF and the FM phases: Measurements sug- 
gest that 7fm w 59 — 62.5 mJ kg~^K~^ is about four to 
six times larger than 7af ~ 10.5 — 16 mJ kg~^K~^>i^ii^ 
Based on these observations, Tu et al. explained the 
transition by a change in entropy of the band electrons 
between the AF and the FM phases. An estimation of 
the free energy shows that these contributions have the 
right order of magnitude to explain the AF-FM transi- 
tion, if one assumes that the electronic densities of states 





FIG. 1: Left: Type II antiferromagnetic ground state struc- 
ture with a nonmagnetic Rh atom at the center and Fe atoms 
with moments of 3.3 /^b at the corner site. Right: Ferromag- 
netic structure with Fe moments of 3.2 at the corners and 
Rh moment of 0.9 ^lb at the center site. 



at the Fermi level do not vary considerably from low tem- 
peratures up to the transition. However, since the FM 
phase cannot be stabilized at low temperatures for stoi- 
chiometric FeRh, 7fm has been measured in slightly more 
iron-rich alloys. With respect to the strong sensitivity of 
the magnetic phase diagram to small departures from 
stoichiometry^ one may doubt that the low tempera- 
ture specific heat can be considered to be independent of 
concentration^ii^ A further drawback considering the 
explanation of Tu et al. arises from the fact that adding 
5 % Iridium boosts 7af by almost an order of magnitude 
to a value of 7af — 101 mJ kg~^K~^>i^*iS The compound 
Fe4g.5Rh45.5Ir5 also undergoes a metamagnetic transition 
with shifted to higher temperatures.-'? Since the rela- 
tion between 7af and 7fm is reversed in this material, 
the previously sketched explanation cannot be applied to 
this case. 

In 1992 Moruzzi and Marcus performed ah initio cal- 
culations using spin polarized density functional theory 
(DFT) in the framework of the local density approxima- 
tion (LDA) and the augmented spherical wave (ASW) 
formalism^iSiSSiSi They calculated the total energy of 
FeRh as a function of volume for different magnetic struc- 
tures. They found that the AF-II spin structure is the 
ground state, whereas the FM structure is another stable 
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Vo AE mpo 

a. u. /atom mRy /^b 
AF FM AF FM 



B N{Ef) 

GPa states/Ry 
AF FM AF FM 



Ref.18 90.2 91.8 1.9 
Ref. ^ 91.9 92.9 2.2 
This work 91.4 93.0 2.5 



2.98 3.15 214 202 37 32 
3.13 3.20 227 244 13 32 
3.18 3.23 197 193 13 29 



TABLE I: Comparison of ab initio results. Here, Vb is the equi- 
hbrium volume of AF and FM phases, AE = Efm — Eaf is 
the energy difference per atom, mpe is the Fe magnetic mo- 
ment, B is the bulk modulus, and N{Ep) the density of states 
per formula unit at the Fermi level. The moment on the Rh 
sites is 1.0 /xb in the FM phase and zero in the AF phase in 
all calculations. 



solution with higher energy and a larger volume. For the 
moments they obtained /ipc = 2.98 /ib, /iRh — 1-02 /ib 
(AF phase) and /xpo = 3.15 /xb (FM phase), in agreement 
with experimental results. 
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In this work, we present further ab initio total energy 
calculations of stoichiometric a-FeRh using the ASW 
formalism and the generalized gradient approximation 
(GGA).^^ The corresponding code also allows to eval- 
uate noncoUinear alignment of spins (see Ref.l23! and ref- 
erences therein) . For many cases LDA gives a reasonable 
description of the ground state properties of solids. How- 
ever, in some cases it predicts the wrong ground state as 
e. g. for iron24. For Fe, the hep structure is found to 
be 10 mRy lower in energy than the bee phase. Because 
FeRh has also a bee ground state lattice structure (un- 
like pure Rh), it is necessary to go beyond LDA by us- 
ing GGA. We have considered different possible magnetic 
ground state structures presented in section ^ Based 
then on the specific energetic order of resulting energy 
versus volume curves, we discuss a new mechanism for 
the temperature driven metamagnetic transition in FeRh. 
In contrast to prior explanations, our model does not 
rely on ground state properties or low temperature data 
alone. Instead, we propose that thermal excitations at 
finite temperatures are the driving force for the transi- 
tion. This is demonstrated in section IIIII on the basis 
of Monte Carlo simulations of an Ising-type spin model, 
showing that a competition between AF Fe-Fe exchange- 
interactions and a nonmagnetic Rh state on the one hand 
and FM Fe-Rh interactions on the other is sufficient to ex- 
plain the metamagnetic transition of FeRh. Inclusion of 
spatial degrees of freedom fsection llV|l by adding pair po- 
tentials and assuming a linear variation of the exchange 
parameter with the lattice parameter proves that this 
model explanation is in accordance with the ab initio 
data and leads furthermore to a nearly quantitative de- 
scription of experimental details. 



FIG. 2: Total energy and magnetic moment versus atomic 
volume as obtained by ab initio calculations (present work) . 



II. AB INITIO TOTAL ENERGY 
CALCULATIONS 

We have performed ab initio calculations by using the 
ASW-method^^ and GGA. Relativistic effects are in- 
cluded in the scalar relativistic approximation. The ba- 
sis wave functions of Fe and Rh atoms include spd and / 
states which are sufficient to obtain the correct magnetic 
behavior of 7-iron. We assume that the AF like spin 
structure can best be described by a spin-spiral (Ref. 
123) with wave vector (0.5,0.5,0.5) in units of 27r/a. As 
first step, we optimized the volume of the AF-II state 
with equal ASA-radii for both types of atoms. For the 
resulting equilibrium volume, the total energy is then 
minimized with respect to the ratio of the ASA-spheres 
leading to rRh/Te = 1-15. In order to investigate the in- 
fluence of the Rh moment in the FM phase we have also 
used the fixed spin moment method^ to restrict mRh to 
zero. Tabled contains the calculated equilibrium proper- 
ties for AF-II and FM phases which are in fair agreement 
with previous calculations. Compared to previous non- 
relativistic ASW LDA calculations (Moruzzi et alm^), we 
obtain equilibrium volumes which are about 1-2 % larger, 
being typical for GGA. In contrast to the cases of pure 
ironSi and FeNi Invar alloys^, the magnetic energy dif- 
ferences are only slightly influenced by the gradient cor- 
rections. Also, no evidence was found that a noncoUinear 
structure could be lower in energy than the previously 
found coUinear AF and FM spin arrangements. With re- 
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spect to the density of states at the Fermi level our results 
compare better with results of Szajek et al2^ than with 
results of Ref. Since the latter results were obtained 
by using LDA ASW (with a spd basis) , we conclude that 
the density of states at the Fermi level is very sensitive to 
computational details. The values obtained for the bulk 
moduli, i?AF and B-pM, are lower compared to those of 
Ref. [l^ (where Baf > -SpM was obtained in contrast to 
Bfm > Saf in Ref.m. 

The calculated total energy curves are shown in Fig.|21 
The large energetic difference between the usual FM 
phase and the hypothetic FM phase with zero rhodium 
moment implies that a finite Rh moment plays an im- 
portant role for the stability of the FM phase. It seems 
unlikely that a magnetic field at the Rh sites induced 
by the surrounding iron atoms is responsible for the ap- 
pearance of a Rh moment in the FM phase, as was pro- 
posed in previous discussions. However, a strong ferro- 
magnetic exchange interaction between the Rh and the 
iron atoms that overrides an antiferromagnetic exchange 
between next nearest neighbor iron sites, would explain 
how the existence of a Rh moment can help to stabilize 
the FM phase. We have evidence as shown in the fol- 
lowing sections, that a competition between a low lying 
nonmagnetic Rh state and another one with higher en- 
ergy and finite moment (which can benefit from exchange 
with ferromagnetic iron neighbors) is the main reason for 
the metamagnetic transition. 
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2.13 


-1.00 






TABLE 11: Parameters for the spin- analogous model (in 
mRy). 

Si — and the magnetic 5^ = ±1 states. For Fe we 
choose a large positive value in order to suppress the 
Si = state. For Rh we choose a negative value leading 
to a nonmagnetic ground state. The second term con- 
tains the exchange parameters Jik which depend only on 
the type of atoms located at the sites i and k. In the 
case of ordered equiatomic FeRh, we have only three dif- 
ferent parameters: Jp"po, JpcRh ^^RhRh- The first one 
is chosen to be negative in order to accomplish an AF 
ground state. The second is taken to be large and pos- 
itive as outlined in the previous section. The third one, 
for the sake of simplicity, is set equal to zero. The choice 
for JpcFc fixed by the Neel-temperature Tn of the AF 
phase assuming that no transition to a FM phase takes 
place. Tn can be determined from the P-T phase diagram 
by extrapolating the transition line between AF and PM 
phases (occurring at pressures larger than 6 GPa) to zero 
pressure. J^^mi ^'^'^ -^Rh have been chosen to yield re- 
alistic values for Tq and Tm, respectively. The values of 
parameters used in the simulations are given in Table ITU 



III. A SPIN-ANALOGOUS MODEL 

In order to confirm our hypothesis we have constructed 
a model being suitable for an examination of the meta- 
magnetic transition at finite temperatures. The simplest 
way to do this is by means of Monte Carlo simulations 
with a localized spin model. To keep the model tractable, 
we neglect spin wave excitations and restrict ourselves to 
Ising spins. This is justified because the transition takes 
place between two ordered structures and both phases 
have collinear spin structures. For a description of the 
nonmagnetic Rh state in addition to the magnetic Fe and 
Rh states, we chose a spin-1 Ising model, where the spin 
variables can take the values Si = — 1,0, +1. The spins 
are located on a bcc lattice with nearest and next near- 
est neighbor interactions. Depending on their positions, 
we distinguish between Fe or Rh sites, where each type 
occupies a simple cubic sublattice, corresponding to an 
ordered equiatomic alloy. The interaction parameters de- 
pend then on the type of sites involved. This situation 
can be described by the following Hamiltonian, 

H = -Y,D,Sl- J^kS^Su. (1) 

i <nn,nnn> 

Without the assumption of different types of atoms, 
Hamiltonian ^ is also known as the Blume-Capel 
modelM^ The first term separates the nonmagnetic 



A. Details of computation 

The evaluation of thermodynamic properties of Q is 
done on the basis of Monte Carlo simulations according 
to the Metropolis schemc"^^ using a sequential update. 
Interesting quantities like magnetization or magnetic mo- 
ment are computed and summed up every 10 to 20 lattice 
sweeps, which ensures that the evaluated lattice config- 
urations are sufficiently uncorrelated. Furthermore, we 
discard the first 20 000 lattice sweeps in order to allow 
the system to reach thermal equilibrium before comput- 
ing averages. In order to speed this up, we have also used 
the final configuration of the last run to initialize the 
simulation for the next temperature. Simulations which 
involve a phase transition are performed twice, with in- 
creasing and decreasing temperatures, in order to assure 
that thermal equilibrium has been reached. 

The computed AF and FM phases are metastable, i. e. 
they are separated by a large energy barrier, which arises 
from the fact that in the transition states a considerable 
amount of FM domains have to be created in the AF 
phase and vice versa. So the standard algorithm is un- 
likely to overcome this barrier and as a result the meta- 
magnetic transition might not be seen at all. Instead 
the phases have to be overheated or undercooled before 
transforming, which results in a large hysteresis or irre- 
versible behavior which makes it difficult to obtain reli- 
able information about the transition point. Therefore, 
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it is necessary to modify the algorithm in order to allow 
for a direct jump to the other phase by circumventing 
the energy barrier with a global update step, where all 
spins will be updated at once. This algorithm ideally con- 
nects equilibrium configurations of the AF phase directly 
with equilibrium configurations of the FM phase while 
ensuring that the entropy difference between the states 
is correctly reproduced. This can be done by choosing a 
unique mapping between lattice configurations of the AF 
and FM phases, respectively. Or, more general, the se- 
lection probability of a specific target configuration must 
be the same as the selection probability of the previous 
start configuration in the backward direction. 

Since it is a priori not clear how equilibrium configu- 
rations will look like at finite temperatures, we simply 
use an update scheme which connects the ground state 
configurations of both phases and thus works at least 
at low temperatures. In the vicinity of the transition 
temperature, this algorithm might not reproduce equi- 
librium states for the trial configuration, because nature 
and amount of excitations are presumably different in 
both phases. Since off-equilibrium means that most of 
the trial states are too high in energy, the probability 
that the trial state is by chance close enough to an equi- 
librium state to be accepted, decreases for larger system 
sizes. 

In order to obtain a trial configuration, we divide the 
system into Rh and Fe sublattices. The Fe sublattice 
is again divided into two sublattices according to spin 
up and spin down positions (see Fig. For each of 
these sublattices, it is decided randomly (with probabil- 
ity 1/2), whether the corresponding sublattice is flipped 
as a whole. For the Rh sublattice another random num- 
ber decides whether the spin values Si £ {—1,0,-1-1} are 
rotated clockwise or counter-clockwise (i.e. —1 becomes 
0, becomes -fl, +1 becomes —1, or vice versa). So, 
each trial configuration is chosen with the same prob- 
ability. Afterwards the energy difference AH between 
the present and the previous configuration is computed 
and the new configuration is accepted with probability 
min(exp(— AiJ/fceT), 1) assuring detailed balance. This 
global update step is performed after each complete lat- 
tice sweep. 

We performed for each temperature between 100 000 
and 1 000 000 (around Tm) lattice sweeps for different lin- 
ear system sizes L = 6 to 16. However, our global update 
scheme does only show a metamagnetic transition within 
reasonable simulation times up to a system size of L = 10. 
Looking for a different way to determine the transition 
point which would also work for larger system sizes, we 
estimate the free energy by integrating the specific heat. 
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^^(r) = E{T) - T 1^5(0) + £^dT^ 



(2) 



In the computation we fit the simulated results for the 
specific heat divided by temperature, C/T, and the in- 
ternal energy E with 10th order polynomials, which can 
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FIG. 3: AF order parameter (diamonds) and magnetization 
(circles) as a function of temperature. The magnetization is 
obtained by multiplying the spin values with a value of 3.0 /ib 
for Fe and 1.0 hb for Rh before averaging. The simulated 
system size is L = 10. 



be integrated analytically. For the system characterized 
by the Hamiltonian (^), we neglect the entropy contri- 
bution at zero temperature >S'(0), since the ground state 
spin structure in both phases is nondegenerate except 
for systems with spin inversion symmetry. Furthermore, 
concerning the estimation of the free energy finite size 
effects are not expected to affect the results, since phase 
transitions are not encountered during these simulations. 
For the rest of the calculations a comparison of the re- 
sults for smaller and larger systems sizes (as can be seen 
in the upcoming figures) reveals that the main issues of 
this paper are also not affected by the restricted system 
sizes. 



B. Computational results 

The order parameter of the AF phase is the staggered 
magnetization for the AF-II spin structure, i. e. the sum 
(of the absolute values) of the staggered magnetizations 
of the two simple cubic sublattices that constitute the bcc 
lattice structure. The staggered magnetization of a sim- 
ple cubic (NaCl) lattice is defined as sum of spins multi- 
plied with a sign which alternates depending on whether 
the corresponding spin occupies a Na or a CI position. 
The order parameter of the FM phase is given by the 
magnetization of the lattice. The variation of both or- 
der parameters with temperature is shown in Fig. |21 for 
a system of size L = 10. At low temperatures the stag- 
gered magnetization approaches a maximum value of 0.5 
due to the fact that the Rh sublattice has no moment 
(Fig. 21) and therefore does not contribute to the sum. 
At Tm = 268 K the staggered magnetization abruptly 
drops to zero, whereas magnetic moment and magneti- 
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FIG. 4: Mean magnetic moment of the Rh atoms. 



zation increase close to their saturation values. Above 
Tc = 610 K ferromagnetism breaks down and the system 
becomes paramagnetic. At the same time, the average 
moment of the Rh atoms falls down to a value of 0.2 /ie- 
This is also the reason, why the decrease in the magneti- 
zation appears unusually sharp. From the data, a phase 
transition of second as well as of first order seems to be 
possible. Accordingly, it is known that for negative val- 
ues of D the Blume-Capel model can show both kinds of 
phase transitions separated by a tricritical point in the 
T-D-planei^SiSi^ Hence in order to safely determine the 
nature of this phase transition further calculations are 
necessary. 

A close look at the specific heat (Fig. O helps to ex- 
plain the occurrence of a metamagnetic transition. At 
first sight it seems as if the specific heats of FM and AF 
phases do not differ very much below 350 K. Above this 
temperature the specific heat of the AF phase increases 
more rapidly with temperature until around 570 K the 
overheated AF phase is not stable anymore in the sim- 
ulations and transforms to the FM phase. However, the 
inset shows that starting around 100 K the specific heat of 
the FM phase is enhanced compared to the AF heat. We 
explain this enhancement by a Schottky-type anomaly 
which adds up to the excitations from spin fiips. Schot- 
tky anomalies are observed in systems with two levels 
separated by a small energy barrier. In fact a crossover 
from magnetic to nonmagnetic Rh atoms is conform with 
this picture. In the FM phase the Rh atoms have a mo- 
ment being ferromagnetically aligned to the Fe moments, 
because the loss of the moment would correspond to an 
energy loss of eight times the exchange constant. This 
amount is diminished by the energy gain due to the Drh 
term, which is smaller than 8 x Jp^^h but larger than 
4 X JpeRh' since the ground state is otherwise not antifer- 
romagnetic. Increasing fluctuations of the magnetization 
of the Fe sublattice cause this energy difference to de- 



FIG. 5: Specific heat as a function of temperature for the FM 
and the AF phases (L — 16) and for simulations on smaller 
lattices using the global MC step (L = 10). Inset: The specific 
heat of the AF and FM phase between and 500 K. Clearly 
visible is a Schottky-type enhancement below 300 K in the FM 
phase, which is responsible for the metamagnetic transition. 
The sharp peak around 270 K corresponds to the metamag- 
netic transition, which is observed in the small system due to 
the global MC step. 
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FIG. 6: Free energy of the FM and AF phases obtained by 
integrating C/T, with the transition point at Tm ~ 266 K. 



crease and finally lead to a breakdown of the average Rh 
moment at Tq- In the AF phase magnetic Rh atoms 
can be excited at the expense of the energy Z^Rh- There 
is, however, no gain in energy due to the exchange in- 
teraction, since the contributions from the AF Fe atoms 
cancel at the Rh site. This corresponds to a much larger 
energy difference between magnetic and nonmagnetic Rh 
states compared to the FM phase and does not lead to 
an appreciable contribution to the specific heat. 
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This view is further supported by a comparison of the 
free energy of both phases shown in Fig. |H| The curves 
intersect at T w 266 K, which is in excellent agreement 
with the simulated = 268 K. This accuracy could be 
achieved, because the data for the specific heat obtained 
from the simulations have only little spread and the fitted 
curves interpolate the data points perfectly. From the 
difference between the internal energies of both phases 
at the transition point, the entropy jump at is deter- 
mined to l^S{T^) = 5.2 JKg-iK-i, which is only 30 to 
45 % of the experimental values (but of the right order 
of magnitude). But, one has to bear in mind that with 
the choice of an Ising model we neglected the possibility 
of non-collinear moments. For the Rh moments this can 
only occur in the FM phase, and would therefore con- 
tribute to AS* (Tin). A second point is that the weight of 
the Si — state is the same as for each magnetic state 
Si — ±1. This is a natural choice for the spin-1 Ising 
model; but for the real system this is somehow arbitrary, 
because we have no information about the electronic ori- 
gin of both Rh states. 



IV. AN EXTENDED MODEL 

Since the simple spin Hamiltonian can reproduce 
a metamagnetic transition as observed in a-FeRh, it re- 
mains an interesting question whether other outstanding 
properties of this alloy, as the large volume increase at 
Tni, can also be explained. Furthermore, it has still to be 
proven whether our spin analogy is a good approximation 
to the ah initio results in the sense that it has compara- 
ble low temperature properties. In order to check this 
we have to extend 1^ for a description of elastic and 
magneto- volume properties: 

= - ^ A - ^ Jikink) Si Sk 

i <nn,nnn> 
+ Yl ^nn(nfc) + Kn„(rifc) . (3) 



For V{rik) we use simple pair potentials of the Lennard- 
Jones type: 



,12 / \ 6 

(J \ la 



Tik 



Tik 



(4) 



Since in general the lattice structure of a Lennard-Jones 
system is closely packed, two different pair potentials for 
nearest and next-nearest neighbors have to be used in 
order to stabilize the bcc structure. The potentials, how- 
ever, do neither distinguish between different atom types 
nor between the different spin states, as has been done 
in previous simulations of related materials like Fe-Ni 
Invar or Y{Mwx^h-x)2^^^ The use of Lennard-Jones 
potentials is far from being optimum for metals, but has 
numerical advantages that enable us to speed up the cal- 
culations substantially. It is then sufhcient to choose the 



(2.6 A) J^S-Us.oA) j;^s;'^,(3.oA) 



>fcBT -11.1 



2.10 



-1.04 







ajg,V/ar 9Jga?J9r dJ^^Jdr 



1.97 



1.58 







25.23 



2.32 



25.23 



2.67 



TABLE III: Magnetic and elastic parameters for the extended 
spin model (energies in mRy, distances in A). 



parameters so that basic elastic properties like low tem- 
perature lattice constant, bulk modulus or thermal ex- 
pansion are reproduced. 

Another change compared to is that the exchange 
parameter is now taken to be a function of interatomic 
distance. For simplicity, we assume a linear distance de- 
pendence: 



Jikink) 



Jii 



dJik 
dr 



Tik 



(5) 



The values for Ji]^ which we take in the extended calcu- 
lations are roughly the same as in Table |n] The relation 
between the derivatives of J^'^Yih "^FcFo '^i^h respect to 
the interatomic distance was determined by the relation 
of the pressure derivatives of the Curie and Neel tem- 
peratures, respectively, which have been obtained from 
the experimental phase diagram by assuming a linear de- 
pendence between dJik/dr and dTc^^/dp. The absolute 
values have been adapted to reproduce the volume jump 
at Tin. An exchange interaction of this form has also 
been used to describe magneto-volume effects in Fe-Ni 
Invari^L^^i^ In this case, the derivative of the exchange 
constant dJik /dr was a factor of 2 to 20 larger than in the 
present work. The values of parameters are summarized 
in Table ITTTl 



A. Details of computation 

For the evaluation of Hamiltonian Q we use a 
textbook isothermal-isobaric Monte Carlo method (e. g. 
Ref. |40|) consisting of alternating spin and position up- 
dates for each atom and a global volume update step af- 
ter finishing each lattice sweep. This algorithm has been 
used by the the authors in previous calculations and is 
explained in detail in the corresponding referencesi^^4^Siii 
In order to simulate the metamagnetic transition, an ad- 
ditional global spin update step has to be introduced. We 
use the same algorithm as described in the last section 
for the spin system in connection with a simultaneous 
volume adaption. Since the latter reduces the acceptance 
probability considerably, the global update scheme has to 
be repeated several thousands of times. It is not practi- 
cable to use a new spin configuration for each trial step, 
because the evaluation of the energy is comparatively 
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FIG. 7: Internal energy and moment as a function of the 
atomic volume obtained for Hamiltonian @ at T = 4K. The 
simulations have been performed for fixed spin structures as 
well as for a freely relaxed spin system. 

time consuming. On the other hand, the new energy 
after solely rescaling the volume can be calculated very 
quickly, since due to the use of Lennard-Jones potentials 
the energy can be written as a function of integral pow- 
ers of the lattice parameter. Therefore, we choose a trial 
spin configuration as described before and compute the 
new energy. Then we attempt a previously fixed number 
N of Metropolis steps, each with a newly chosen volume. 
If one step is accepted, we continue with the original spin 
configuration as trial system (and so on) until N steps 
have been made. Since the number of trial steps N has 
been previously fixed, detailed balance is still valid. 

We have performed simulations with system sizes of 
i = 6 to 12. A direct metamagnctic transition could 
only be seen for system sizes up to L = 8. The simulation 
time ranged from 120 000 up to 1000 000 lattice sweeps 
around with values from 1000 to 10 000 for N. As 
before we estimated Gibbs' free energy for zero pressure 
by integrating the specific heat according to Eq. 



B. Computational results 

For a comparison of our model properties with the re- 
sults of ab initio calculations, we calculated in isochoric 
simulations the energy as a function of the volume at low 
temperatures for different fixed spin structures (Fig. [7)). 
Here, we cooled a system of size L — 16 exponentially 
down from 100 K down to 4K. We find good agreement 
with the results of Fig. |21 in the sense that the order of 
magnetic phases is similar. Although the ab initio to- 



FIG. 8: AF order parameter (diamonds) and magnetization 
(circles) as a function of temperature as obtained for the ex- 
tended model Hamiltonian The system size is L = 8. 
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FIG. 9: Mean magnetic moment of the Rh atoms for the 
extended model {L = 8). 



tal energy differences are somewhat larger, this indicates 
that our model Hamiltonian is a qualitatively correct 
description of the mechanisms leading to a metamagnctic 
transition in FcRh. 

As in the simple model without volume-dependent 
terms, we find an abrupt increase of the magnetiza- 
tion in combination with a discontinuous decrease for 
the staggered magnetization with increasing tempera- 
ture (Fig. IS)). Consequently the mean moment at the 
Rh sites (Fig. O also raises sharply around the metam- 
agnctic transition temperature of Tm = 322 K. Around 
the Curie temperature T = 720 K the magnetization de- 
creases more smoothly than for the simple model. This 
may be due to the enhancement of the effective exchange 
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FIG. 10: Volume per atom plotted over temperature as ob- 
tained for the AF and FM phases (L — 12) and by our simu- 
lations with the global update scheme (L = 8). 
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FIG. 11: Adiabatic bulk modulus as a function of temperature 
for the AF and FM phases {L — 12) and with the global 
update scheme (L — 8). The latter data, which show the 
metamagnetic transition, reveal a very large softening of the 
lattice around Tm. 



parameter given by Eq. (jSJ which is caused by the lattice 
expansion. 

As expected from the low temperature calculations, 
the volume of the AF phase is smaller than the volume 
of the FM phase throughout the stability range fFig. llO|l . 
For the freely fluctuating system a volume jump of 0.8 % 
occurs at Tm. Below Tc the volume expansion is reduced 
in the FM phase, which is in qualitative agreement with 
experiment!***. 

The bulk modulus i?s of the AF phase is about 6 % 
larger than the FM bulk modulus throughout the stabil- 



FIG. 12: Specific heat for the extended model as a function 
of temperature for FM and the AF phases (L = 12) and 
for simulations on smaller lattices using the global MC step 
(L = 8). The relationship between the specific heats of AF 
and FM phases below 400 K is shown in the inset. 



ity range fFig. Ill|l . Around Tm and Tc we find a con- 
siderable weakening of the material which is caused by 
the magneto- volume anomalies. The values for Bg at low 
temperatures can be estimated more accurately by using 
a polynomial fit to the E-V curves in Fig. [7| We find 
Bs = 231 GPa for the AF phase and Bs = 217 GPa for 
the FM phase. These values are in good agreement with 
results of isobaric calculations shown in Fig. ^] Com- 
pared with the results of ab initio calculations (Table P) 
the absolute values are too large, whereas these calcu- 
lations do not give a unanimous prediction for the sign 
and the magnitude of the difference A_Bg (with respect to 
the magnetic structures). Experimental measurements^ 
of Young modulus suggest that the bulk modulus of the 
AF phase should in fact be lower than in the FM phase. 
The modulus of the AF phase, however, has only been 
estimated in the vicinity of Tm, where a weakening of the 
material (as a precursor of the transition) is present (as 
in our simulations). 

The specific heats obtained for the extended model 
(Fig. I12|) resemble the findings for the spin-only Hamil- 
tonian (Q, except that now additional contributions due 
to atomic displacements are included. This shows that 
the proposed explanation of the metamagnetic transition 
is still valid. Since the Hamiltonian (|2Jl does not contain 
any kinetic terms, the low temperature value of the spe- 
cific heat is only half of the Dulong-Pctit limit of 3 R, 
where R is the kinetic gas constant. 



C. Estimation of Gibbs free energy 

The calculation of Gibbs free energy G for the model 
with classical motion (O is more complicated than in the 
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FIG. 13: Differences of internal energies and Gibbs free en- 
ergies for FM and AF pliases obtained by integrating Cp/T. 
Tfie transition point is Tm = 324 K. Tlie experimental data 
have been obtained by integrating over the hysteresis loop for 
different temperatures^ 



spin-only case (^3) • First of all absolute values of G can- 
not be given, since the specific heat at zero temperature 
is finite and the resulting entropy would diverge. But 
since the limit for T^O of Cp is 1.5 i? and hence the 
same for all magnetic structures, differences of free ener- 
gies can be computed. In contrast to the previous sec- 
tion, the entropy contribution at zero temperature A5'(0) 
must be considered. As before a contribution from the 
magnetic system can be neglected, while we have to ac- 
count for the differences in the elastic properties of FM 
and the AF phases. Since we only need the entropy at 
zero temperature, we choose an ensemble of harmonic 
oscillators as an approximation for our spatial degrees of 
freedom, neglecting the anharmonicity of the potentials 
and the coupling of the oscillators. Differentiating the 
free energy obtained from the logarithm of the partition 
function leads then to a simple expression for the entropy 
difference. 



AS- = -NkB In (fcAF/fcpM) 



(6) 



where ^af and fcpM are the force constants of the har- 
monic potentials which are estimated from the curva- 
ture of the ground state energy (versus lattice parame- 
ter) curve. For the entropy difference we obtain then: 
AS'(O) = 9.17Jkg"^K-^ Taking this into account, we 
achieve again a rather good value for Tm. From the dif- 
ferences of internal energies at Tm we obtain for the en- 
tropy jump at the transition IS.S{T^) — 15.9 J kg^^K^^, 
which is within the range of experimental results. Com- 
parison of the calculated free energy with experimental 
values for AG shows excellent agreement. The experi- 
mental values have been obtained by a graphical inte- 
gration of the magnetic field expressed as a function of 



the measured magnetization:-- Extrapolation of experi- 
mental values to zero temperature shows that the energy 
difference Ai?(0) between the AF and FM phases is much 
better described by our model parameters than by ah ini- 
tio results, since the latter show that Ai? is one order of 
magnitude too large. This discrepancy has already been 
noticed by Moruzzi and Marcus^^ who relate this to the 
omission of zero point energy corrections in their total 
energy calculations. 



V. CONCLUSIONS 

We propose on the basis of new ah initio results a mech- 
anism for the metamagnetic transition in FeRh at finite 
temperatures. In contrast to previous explanations, our 
model does not rely on a large difference between the 
low temperature specific heat constants of both phases. 
These are expected to be sensitive to external influences 
as both experimental measurements and band structure 
calculations suggest, so that it seems implausible that 
a constant contribution of the given magnitude might 
survive up to room temperature. Instead, we propose 
that the existence of two magnetic states of Rh atoms 
connected with competing FM Fe-Rh and AF Fe-Fe ex- 
change interactions are at the origin of the metamagnetic 
transition. The magneto-volume effects can simply be 
explained on the basis of distance dependent exchange 
parameters. The applicability of this mechanism to the 
Fe-Rh problem has been verified by Monte Carlo model 
calculations, showing that a metamagnetic transition of 
the desired kind does in fact occur and, by extending the 
model, magneto-volume effects and other experimental 
properties can be sufficiently well described. 

As we have pointed out, our explanation is in agree- 
ment with existing experimental data. However, a fur- 
ther check of our model would be a comparison with 
(non-existing) specific heat data from above Tq down to 
very low temperatures for both, the AF and FM phases. 
From this one could then estimate the magnetic contri- 
bution by subtracting the lattice part within the Debye 
approximation, the electronic part and the contribution 
by the anharmonicity of the potentials42il4 So far, the 
specific heat has only been determined in the range from 
100 K to 500 K and only for the nearly stoichiometric 
Rh-rich alloy4S^ Systematic measurements on the Fe- 
rich side with a FM ground state and measurements un- 
der pressure, suppressing the metamagnetic transition 
would yield information whether a Schottky-type exci- 
tation plays an important role, which should show up 
around 200 K in the magnetic contribution to the spe- 
cific heat of ferromagnetic samples. 

Monte Carlo simulations with applied pressure are left 
for future work, since with increasing pressure and hence 
increasing Tm a reliable estimation of the transition tem- 
perature is rather difficult and requires an improvement 
of the global MC step. First tests, however, showed that 
the location of the phase boundaries under pressure is in 
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sufficient agreement with experimental data for pressures 
below 20kbar and above 40kbar, where the metamag- 
netic transition is completely suppressed. The tricritical 
point, if one exists, should be located somewhere between 
20 kbar and 40 kbar for the parameters used here. 
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